Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense.

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Genome-wide investigation of SQUAMOSA promoter binding protein-like genes in Liriodendron and functional characterization of LcSPL2.

The plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play a pivotal role in various developmental processes, including leaf morphogenesis and vegetative to reproductive phase transition. Liriodendron chinense and Liriodendron tulipifera are widely used in landscaping due to their tulip-like flowers and peculiar leaves. However, the SPL gene family in Liriodendron has not been identified and systematically characterized. We systematically identified and characterized the SPL family members in Liriodendron, including phylogeny, gene structure and syntenic analyses. Subsequently, we quantified the expression patterns of LcSPLs across various tissue sites through transcription-quantitative polymerase chain reaction (RT-qPCR) assays and identified the target gene, LcSPL2. Finally, we characterized the functions of LcSPL2 via ectopic transformation. Altogether, 17 LcSPL and 18 LtSPL genes were genome-widely identified in L. chinense and L. tulipifera, respectively. All the 35 SPLs were grouped into 9 clades. Both species had three SPL gene pairs arising from segmental duplication events, and the LcSPLs displayed high collinearity with the L. tulipifera genome. RT-qPCR assays showed that SPL genes were differentially expressed in different tissues, especially. Because LcSPL2 is highly expressed in pistils and leaves, it was selected to describe the SPL gene family of L. chinense by ectopic expression. We showed that overexpression of LcSPL2 in Arabidopsis thaliana resulted in earlier flowering and fewer rosette leaves. Moreover, we observed that overexpression of LcSPL2 in A. thaliana up-regulated the expression levels of four genes related to flower development. This study identified SPL genes in Liriodendron and characterized the function of LcSPL2 in advancing flower development.

The identification and functional characterization of the LcMCT gene from Liriodendron chinense reveals its potenatial role in carotenoids biosyanthesis.

Terpenoids are not only important component of plant floral scent, but also indispensable elements in the formation of floral color. The petals of Liriodendron chinense are rich in tetraterpene carotenoids and release large amounts of volatile monoterpene and sesquiterpene compounds during full blooming stage. However, the mechanism of terpenoid synthesis is not clear in L. chinense. In this study, we identified a LcMCT gene and characterized its potential function in carotenoids biosynthesis. A total of 2947 up-regulated differentially expressed genes (DEGs) were discerned from the transcriptomic data of L. chinense petals, with a significant enrichment of DEGs related to plant hormone signal transduction and terpenoid backbone biosynthesis. After comprehensive analysis on these DEGs, the LcMCT gene was selected for subsequent function characterization. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) results showed that LcMCT was expressed at the highest level in the petals during full blooming stage, suggesting a possible role in carotenoids biosynthesis and volatile terpenoid biosynthesis. Subcellular localization showed that the LcMCT protein was localized in the chloroplast. Overexpression of LcMCT in Arabidopsis thaliana affected the expression levels of MEP pathway genes. Moreover, the MCT enzyme activity and carotenoids contents in transgenic A. thaliana were increased by 69.27% and 15.57%, respectively. These results suggest that LcMCT promotes the biosynthesis of terpenoid precursors via the MEP pathway. Our work lays a foundation for exploring the mechanism of terpenoid synthesis in L. chinense.

Genome-wide identification of GROWTH-REGULATING FACTORs in Liriodendron chinense and functional characterization of LcGRF2 in leaf size regulation.

GROWTH-REGULATING FACTORs (GRFs) play a pivotal role in the regulation of leaf size in plants and have been widely reported in plants. However, their specific functions in leaf size regulation in Liriodendron chinense remains unclear. Therefore, in this study, we identified GRF genes on a genome-wide scale in L. chinense to characterize the roles of LcGRFs in regulating leaf size. A total of nine LcGRF genes were identified, and these genes exhibited weak expression in mature leaves but strong expression in shoot apex. Notably, LcGRF2 exhibited the highest expression level in the shoot apex of L. chinense. Further RT-qPCR assay revealed that the expression level of LcGRF2 gradually decreased along with the leaf development process, and also displayed a gradient along the leaf proximo-distal and medio-lateral axes. Furthermore, overexpression of LcGRF2 in Arabidopsis thaliana resulted in increased leaf size, and significantly up-regulated the expression of genes involved in cell division like AtCYCD3;1, AtKNOLLE, and AtCYCB1;1, indicating that LcGRF2 may influence leaf size by promoting cell proliferation. This work contributes to a better understanding of the roles and molecular mechanisms of LcGRFs in the regulation of leaf size in L. chinense.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories.

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Overexpression of LtuHB6 from Liriodendron tulipifera causes lobed-leaf formation in Arabidopsis thaliana.

Liriodendron tulipifera L. is an ornamental tree species with extraordinarily lobed leaves. However, the mechanisms underlying lobed leaf formation in plants remain unclear. The transcription factor, ARABIDOPSIS THALIANA HOMEBOX 6 (HB6), plays a role in regulating leaf margin development. HB6 is involved in cell division and differentiation of developmental organs and negatively regulates abscisic acid (ABA) signal transmission under external abiotic stress; it is unclear whether HB6 performs a pivotal role in leaf morphogenesis in L. tulipifera. In this study, full-length LtuHB6 from L. tulipifera was heterologously expressed in tobacco and Arabidopsis thaliana; its expression pattern was analyzed to determine its potential role in leaf development. In addition, LtuHB6 is localized in the nucleus and cell membrane of tobacco leaves. The expression of LtuHB6 was highest in mature leaves compared to the other stages of leaf development (bud growth, young leaves, and leaf senescence). Transgenic A. thaliana plants overexpressing LtuHB6 exhibited an abnormal phenotype with lobed leaves. Moreover, LtuHB6 overexpression significantly affected the expression of seven genes related to leaf serration in the initial stage of leaf primordia and altered the expression levels of hormonal genes. Our findings indicate that LtuHB6 is an essential regulatory factor in L. tulipifera lobed-leaf formation and is involved in regulating and responding to hormones. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01254-9.

Genome-wide identification of XTH genes in Liriodendron chinense and functional characterization of LcXTH21.

Liriodendron chinense is a relic tree species of the family Magnoliaceae with multiple uses in timber production, landscape decoration, and afforestation. L. chinense often experiences drought stress in arid areas. However, the molecular basis underlying the drought response of L. chinense remains unclear. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in drought stress resistance. Hereby, to explore the drought resistance mechanism of L. chinense, we identify XTH genes on a genome-wide scale in L. chinense. A total of 27 XTH genes were identified in L. chinense, and these genes were classified into three subfamilies. Drought treatment and RT-qPCR analysis revealed that six LcXTH genes significantly responded to drought stress, especially LcXTH21. Hence, we cloned the LcXTH21 gene and overexpressed it in tobacco via gene transfer to analyze its function. The roots of transgenic plants were more developed than those of wild-type plants under different polyethylene glycol (PEG) concentration, and further RT-qPCR analysis showed that LcXTH21 highly expressed in root compared to aboveground organs, indicating that LcXTH21 may play a role in drought resistance through promoting root development. The results of this study provide new insights into the roles of LcXTH genes in the drought stress response. Our findings will also aid future studies of the molecular mechanisms by which LcXTH genes contribute to the drought response.

Liriodendron chinense LcMAX1 regulates primary root growth and shoot branching in Arabidopsis thaliana.

Strigolactones (SLs) play prominent roles in regulating shoot branching and root architecture in model plants. However, their roles in non-model (particularly woody) plants remain unclear. Liriodendron chinense is a timber tree species widely planted in southern China. The outturn percentage and wood quality of L. chinense are greatly affected by the branching characteristics of its shoot, and the rooting ability of the cuttings is key for its vegetative propagation. Here, we isolated and analyzed the function of the MORE AXILLARY GROWTH 1 (LcMAX1) gene, which is involved in L. chinense SL biosynthesis. RT-qPCR showed that LcMAX1 was highly expressed in the roots and axillary buds. LcMAX1 was located in the endoplasmic reticulum (ER) and nucleus. LcMAX1 ectopic expression promoted primary root growth, whereas there were no phenotypic differences in shoot branching between transgenic and wild-type (WT) A. thaliana plants. LcMAX1 overexpression in the max1 mutant restored them to the WT A. thaliana phenotypes. Additionally, AtPIN1, AtPIN2, and AtBRC1 expressions were significantly upregulated in transgenic A. thaliana and the max1 mutant. It was therefore speculated that LcMAX1 promotes primary root growth by regulating expression of auxin transport-related genes in A. thaliana, and LcMAX1 inhibits shoot branching by upregulating expression of AtBRC1 in the max1 mutant. Altogether, these results demonstrated that the root development and shoot branching functions of LcMAX1 were similar to those of AtMAX1. Our findings provide a foundation for obtaining further insights into root and branch development in L. chinense.

Predicting the impact of climate change on the distribution of two relict Liriodendron species by coupling the MaxEnt model and actual physiological indicators in relation to stress tolerance.

Climate change has a crucial impact on the distributions of plants, especially relict species. Hence, predicting the potential impact of climate change on the distributions of relict plants is critical for their future conservation. Liriodendron plants are relict trees, and only two natural species have survived: L. chinense and L. tulipifera. However, the extent of the impact of future climate change on the distributions of these two Liriodendron species remains unclear. Therefore, we predicted the suitable habitat distributions of two Liriodendron species under present and future climate scenarios using MaxEnt modeling. The results showed that the area of suitable habitats for two Liriodendron species would significantly decrease. However, the two relict species presented different habitat shift patterns, with a local contraction of suitable habitat for L. chinense and a northward shift in suitable habitat for L. tulipifera, indicating that changes in environmental factors will affect the distributions of these species. Among the environmental factors assessed, May precipitation induced the largest impact on the L. chinense distribution, while L. tulipifera was significantly affected by precipitation in the driest quarter. Furthermore, to explore the relationship between habitat suitability and Liriodendron stress tolerance, we analyzed six physiological indicators of stress tolerance by sampling twelve provenances of L. chinense and five provenances of L. tulipifera. The composite index of six physiological indicators was significantly negatively correlated with the habitat suitability of the species. The stress tolerance of Liriodendron plants in highly suitable areas was lower than that in areas with moderate or low suitability. Overall, these findings improve our understanding of the ecological impacts of climate change, informing future conservation efforts for Liriodendron species.

Stable reference gene selection for quantitative real-time PCR normalization in passion fruit (Passiflora edulis Sims.).

BACKGROUND: Passiflora edulis is a tropical fruit with high nutrient and medicinal values that is widely planted in southern China. However, the molecular biology of P. edulis has not been well studied. There are few reports regarding the choice of reference genes for gene expression studies of passion fruit. METHODS AND RESULTS: By using three algorithms, implemented in geNorm, NormFinder and BestKeeper, we have selected ten candidate reference genes to explore their transcriptional expression stability in various tissues and under cold stress conditions. EF1 and HIS were stably expressed in five tissues. Ts and OTU were stably in vegetative organs. 50 S and Liom were stably in reproductive organs. The transcriptional abundance of EF1 and UBQ was stable in cold-treated and recovery treated leaf samples of P. edulis. In all samples, EF1 and Ts exhibited the highest expression stability. Evaluation of selected genes using simple statistical methods (ANOVA and post hoc analysis). Overall, EF1 emerged as the optimum reference gene for qRT-PCR normalize in P. edulis. In addition, the qRT-PCR analysis revealed that expression of ICE1 increases with the duration of cold treatment. CONCLUSIONS: In this study, we successfully screened stable reference genes from 10 candidates in P. edulis and verified the results by analyzing the expression level of ICE1. The results provide reliable and effective reference genes for future research on gene expression analysis in P. edulis, and lay a foundation for follow-up research on functional genes in P. edulis.

The Roles of microRNA-Long Non-coding RNA-mRNA Networks in the Regulation of Leaf and Flower Development in Liriodendron chinense.

The leaf and the flower are vital plant organs owing to their roles in photosynthesis and reproduction. Long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and transcription factors (TFs) are very important to the development of these organs. Liriodendron chinense is a common ornamental tree species in southern China with an unusual leaf shape and tulip-like flowers. The genetic mechanisms underlying leaf and flower development in L. chinense and the miRNA-lncRNA-TF regulatory networks are poorly studied. Through the integration and analysis of different types of sequencing data, we identified the miRNA-lncRNA-TF regulatory networks that were related to leaf and flower development. These networks contained 105 miRNAs, 258 lncRNAs, 393 TFs, and 22 endogenous target mimics. Notably, lch-lnc7374-miR156h-SPL3 and lch-lnc7374-miR156j-SPL9 were potential regulators of stamen and pistil development in L. chinense, respectively. miRNA-lncRNA-mRNA regulatory networks were shown to impact anther development, male and female fertility, and petal color by regulating the biosynthesis of phenylpropanoid metabolites. Phenylpropanoid metabolite biosynthesis genes and TFs that were targeted by miRNAs and lncRNAs were differentially expressed in the leaf and flower. Moreover, RT-qPCR analysis confirmed 22 differentially expressed miRNAs, among which most of them showed obvious leaf or flower specificity; miR157a-SPL and miR160a-ARF module were verified by using RLM-RACE, and these two modules were related to leaf and flower development. These findings provide insight into the roles of miRNA-lncRNA-mRNA regulatory networks in organ development and function in L. chinense, and will facilitate further investigation into the regulatory mechanisms of leaf and flower development in L. chinense.

Overexpression of Liriodendron tulipifera JAG Gene (LtuJAG) Changes Leaf Shapes in Transgenic Arabidopsis thaliana.

In Arabidopsis thaliana, JAGGED (JAG) is a transcription inhibitor that controls the development of leaf polarity and regulates the expression of genes controlling lateral organ formation. Liriodendron tulipifera is an ornamental tree with extraordinary tulip-shaped flowers and goose web-like leaves, this is one of the suitable plants for morphological development research. To investigate the potential functions of the LtuJAG gene, we isolated the full-length LtuJAG from L. tulipifera, transferred it into A. thaliana via agrobacterium-mediated transformation, and monitored its expression pattern. Subcellular localization showed that LtuJAG was located in the nucleus. RT-qPCR assays indicated that LtuJAG was expressed mainly in leaf buds and flowers, but not in mature leaves and stems. GUS staining results showed that LtuJAG was expressed in the shoot apical meristem (SAM). Overexpressing LtuJAG changed A. thaliana leaf shapes, causing a moderate serration and a slight asymmetric distribution in the medio-lateral and proximal-distal axes. Ectopic expression of LtuJAG induced the expression of lateral organ boundary suppressors JAGGED LATERAL ORGANS (JLO) and ARABIDOPSIS THALIANA HOMEOBOX1 (ATH1). It also repressed the expression of the apical meristem suppressor class-1 KNOX gene (KNOX I) and altered endogenous hormone levels. Our results suggest that LtuJAG plays a role in negatively regulating leaf polarity formation in L. tulipifera.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments.

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

Genome-wide identification and characterization of LcCCR13 reveals its potential role in lignin biosynthesis in Liriodendron chinense.

Introduction: Wood formation is closely related to lignin biosynthesis. Cinnamoyl-CoA reductase (CCR) catalyzes the conversion of cinnamoyl-CoA to cinnamaldehydes, which is the initiation of the lignin biosynthesis pathway and a crucial point in the manipulation of associated traits. Liriodendron chinense is an economically significant timber tree. Nevertheless, the underlying mechanism of wood formation in it remains unknown; even the number of LcCCR family members in this species is unclear. Materials and Results: This study aimed to perform a genome-wide identification of genes(s) involved in lignin biosynthesis in L. chinense via RT-qPCR assays and functional verification. Altogether, 13 LcCCR genes were identified that were divided into four major groups based on structural and phylogenetic features. The gene structures and motif compositions were strongly conserved between members of the same groups. Subsequently, the expression patterns analysis based on RNA-seq data indicated that LcCCR5/7/10/12/13 had high expression in the developing xylem at the stem (DXS). Furthermore, the RT-qPCR assays showed that LcCCR13 had the highest expression in the stem as compared to other tissues. Moreover, the overexpression of the LcCCR13 in transgenic tobacco plants caused an improvement in the CCR activity and lignin content, indicating that it plays a key role in lignin biosynthesis in the stems. Discussion: Our research lays a foundation for deeper investigation of the lignin synthesis and uncovers the genetic basis of wood formation in L. chinense.

Genome-wide survey and identification of AP2/ERF genes involved in shoot and leaf development in Liriodendron chinense.

BACKGROUND: Liriodendron chinense is a distinctive ornamental tree species due to its unique leaves and tulip-like flowers. The discovery of genes involved in leaf development and morphogenesis is critical for uncovering the underlying genetic basis of these traits. Genes in the AP2/ERF family are recognized as plant-specific transcription factors that contribute to plant growth, hormone-induced development, ethylene response factors, and stress responses. RESULTS: In this study, we identified 104 putative AP2/ERF genes in the recently released L. chinense genome and transcriptome database. In addition, all 104 genes were grouped into four subfamilies, the AP2, ERF, RAV, and Soloist subfamilies. This classification was further supported by the results of gene structure and conserved motif analyses. Intriguingly, after application of a series test of cluster analysis, three AP2 genes, LcERF 94, LcERF 96, and LcERF 98, were identified as tissue-specific in buds based on the expression profiles of various tissues. These results were further validated via RT-qPCR assays and were highly consistent with the STC analysis. We further investigated the dynamic changes of immature leaves by dissecting fresh shoots into seven discontinuous periods, which were empirically identified as shoot apical meristem (SAM), leaf primordia and tender leaf developmental stages according to the anatomic structure. Subsequently, these three candidates were highly expressed in SAM and leaf primordia but rarely in tender leaves, indicating that they were mainly involved in early leaf development and morphogenesis. Moreover, these three genes displayed nuclear subcellular localizations through the transient transformation of tobacco epidermal cells. CONCLUSIONS: Overall, we identified 104 AP2/ERF family members at the genome-wide level and discerned three candidate genes that might participate in the development and morphogenesis of the leaf primordium in L. chinense.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense.

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

Isolation, expression, and functional analysis of the geranylgeranyl pyrophosphate synthase (GGPPS) gene from Liriodendron tulipifera.

Terpenoids are important secondary metabolites in plants and are involved in stress responses and pollinator attraction. Geranylgeranyl pyrophosphate synthase (GGPPS) is a key synthase in the 2C-methyl-D-erythritol-4-phosphate (MEP) pathway of terpenoid synthesis, catalyzing the synthesis of diterpenoids. Liriodendron tulipifera is a nectar plant in North America. Little is known about the key genes involved in the biosynthetic pathways of terpenoids, the precursors of most compounds related to nectar, fragrance and coloring in flowers in L. tulipifera. In this study, the LtuGGPPS2 gene and its promoter (LtuGGPPS2-pro) were cloned from L. tulipifera. The results of sequence alignment showed that the LtuGGPPS2 gene is highly homologous to GGPPS genes of other plants. Subcellular localization analysis showed that the LtuGGPPS2 protein localizes to chloroplasts, suggesting that the LtuGGPPS2 gene is probably related to carotenoid and chlorophyll synthesis. Based on tissue expression profiles revealed by RT-qPCR, the expression level of the LtuGGPPS2 gene was highest in petals. These results were consistent with the changes in volatile and nonvolatile terpenoids in the flowers of L. tulipifera. GUS staining to examine the LtuGGPPS2 promoter indicated that it is responsive to hormones. Overexpression of the LtuGGPPS2 gene increased the carotenoid content and GGPPS enzyme activity in Arabidopsis thaliana, indicating that LtuGGPPS2 is the key terpenoid synthase in the flowers of L. tulipifera. Our findings lay a foundation for further functional analysis of the LtuGGPPS2 gene and deeper investigation of the terpenoid biosynthetic pathway in L. tulipifera.

Genome-wide identification of MIKC-type genes related to stamen and gynoecium development in Liriodendron.

The organogenesis and development of reproductive organs, i.e., stamen and gynoecium, are important floral characteristics that are closely related to pollinators and reproductive fitness. As a genus from Magnoliaceae, Liriodendron has only two relict species: L. chinense and L. tulipifera. Despite the similar flower shapes of these species, their natural seed-setting rates differ significantly, implying interspecies difference in floral organogenesis and development. MADS-box genes, which participate in floral organogenesis and development, remain unexplored in Liriodendron. Here, to explore the interspecies difference in floral organogenesis and development and identify MADS-box genes in Liriodendron, we examined the stamen and gynoecium primordia of the two Liriodendron species by scanning electron microscopy combined with paraffin sectioning, and then collected two types of primordia for RNA-seq. A total of 12 libraries were constructed and 42,268 genes were identified, including 35,269 reference genes and 6,999 new genes. Monoterpenoid biosynthesis was enriched in L. tulipifera. Genome-wide analysis of 32 MADS-box genes was conducted, including phylogenetic trees, exon/intron structures, and conserved motif distributions. Twenty-six genes were anchored on 17 scaffolds, and six new genes had no location information. The expression profiles of MIKC-type genes via RT-qPCR acrossing six stamen and gynoecium developmental stages indicates that the PI-like, AG/STK-like, SEP-like, and SVP-like genes may contribute to the species-specific differentiation of the organogenesis and development of reproductive organs in Liriodendron. Our findings laid the groundwork for the future exploration of the mechanism underlying on the interspecific differences in reproductive organ development and fitness in Liriodendron.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense.

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Cloning, Characterization and Functional Analysis of the LtuPTOX Gene, a Homologue of Arabidopsis thalianaIMMUTANS Derived from Liriodendron tulipifera.

Flower colour and colour patterns are crucial traits for ornamental species; thus, a comprehensive understanding of their genetic basis is extremely significant for plant breeders. The tulip tree (Liriodendron tulipifera Linn.) is well known for its flowers, odd leave shape and tree form. However, the genetic basis of its colour inheritance remains unknown. In this study, a putative plastid terminal oxidase gene (LtuPTOX) was identified from L. tulipifera based on multiple databases of differentially expressed genes at various developmental stages. Then, the full-length cDNA of LtuPTOX was derived from tepals and leaves using RACE (rapid amplification of cDNA ends) approaches. Furthermore, gene structure and phylogenetic analyses of PTOX as well as AOXs (alternative oxidases), another highly similar homologue in the AOX family, were used to distinguish between the two subfamilies of genes. In addition, transient transformation and qPCR methods were used to determine the subcellular localization and tissue expression pattern of the LtuPTOX gene. Moreover, the expression of LtuPTOX as well as pigment contents was investigated to illustrate the function of this gene during the formation of orange bands on petals. The results showed that the LtuPTOX gene encodes a 358-aa protein that contains a complete AOX domain (PF01786). Accordingly, the LiriodendronPTOX and AOX genes were identified as only paralogs since they were rather similar in sequence. LtuPTOX showed chloroplast localization and was expressed in coloured organs such as petals and leaves. Additionally, an increasing pattern of LtuPTOX transcripts leads to carotenoid accumulation on the orange-band during flower bud development. Taken together, our results suggest that LtuPTOX is involved in petal carotenoid metabolism and orange band formation in L. tulipifera. The identification of this potentially involved gene will lay a foundation for further uncovering the genetic basis of flower colour in L. tulipifera.